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Abstract. We investigate steady states of a quasilinear first order 
hyperbolic partial intcgro-differential equation. The model describes 
the evolution of a hierarchical structured population with distributed 
states at birth. Hierarchical size-structured models describe the dy- 
namics of populations when individuals experience size-specific envi- 
ronment. This is the case for example in a population where individu- 
als exhibit cannibalistic behavior and the chance to become prey (or 
to attack) depends on the individual's size. The other distinctive fea- 
ture of the model is that individuals are recruited into the population 
at arbitrary size. This amounts to an infinite rank integral operator 
describing the recruitment process. First we establish conditions for 
the existence of a positive steady state of the model. Our method 
uses a fixed point result of nonlinear maps in conical shells of Banach 
spaces. Then we study stability properties of steady states for the 
special case of a separable growth rate using results from the theory 
of positive operators on Banach lattices. 

1. Introduction. Classic population models often assume that individu- 
als experience scramble competition. This means that all individuals in the 
population have equal chances in the competition for resources such as food, 
light, space etc., see e.g. [9, 18, 23,[24l[27]. In many species, however, compe- 
tition for resources that determine individual mortality and fertility is based 
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on some hierarchy in the population which is related to individuals size or 
to any other variable that characterizes physiological structure. Significant 
amount of interest has been paid to understand the dynamics of popula- 
tions that exhibit contest competition. Structured population models are a 
useful tool to study intra-specific contest competition. Both (time-)discretc 
(see e.g. H3I2T]) and continuous (see e.g. [TJ |2J M El US ) models have 
been formulated and analysed to this end. Most of the nonlinear models 
in the literature incorporate environmental feedback through some form of 
density dependence in the vital rates. Our goal in this paper is to carry 
out a qualitative analysis of a continuous model which incorporates quite 
general nonlincaritics. 

We study the following quasilincar partial intcgro-diffcrcntial equation 

d d 

-g t P(s,t) + — { 7 {s,P(t))p(s,t)) = -fx(s,E(s,p))p(s,t) 



+ / p(s,y,E(y,p))p(y,t)dy, (1.1) 
Jo 

7 (o,p(i)MM) = o, (i.2) 

E(s,p) = a / w(r)p(r, t) dr + / w(r)p(r,t) dr, (1-3) 



P(t)= / K (s)p(s,t)ds, (1.4) 



with the initial condition 

p(s,0) =p (s). 



This model describes the long term dynamics of a population of a suffi- 
ciently large size living in a closed habitat. The function p = p(s, t) denotes 
the density of individuals of size (or other developmental stage) s at time t. 
It is assumed that individuals may have different sizes at birth and there- 
fore j3(s, y, ■ ) denotes the rate at which individuals of size y "produce" 
individuals of size s. Hence the non-local integral term in Equation p. II) 
represents the recruitment of individuals into the population. 7 denotes the 
size-specific growth rate, while fi stands for the mortality rate. We assume 
that individual growth is also regulated by a weighted population size P, 
for example due to competition. Mortality however, depends on the size- 
specific environment E, for example due to cannibalism. The parameter 
a is related to the strength of the hierarchy in the population. We note 
that if a — 1 our model reduces to the scramble competition model that in 
the special case k = w = 1 was considered in [12] ■ We make the following 
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general assumptions on the model ingredients 

fx E C^m] x [0,oo)), /3 £ C^O.m] x [0,m] x [0,oo)), 

7 £ C 2 ([0, m] x [0, oo)), w, k £ L°°(0, m), (1.5) 

/?, a, iu, > 0, 7, k > 0. 

Notice that we make no requirement that P(s,y, E(y,p)) = if y < s 
although this seems natural from a biological point of view. In the remark- 
able paper [7J , Calsina and Saldaha studied well-posedness of a very general 
size-structured model with distributed states at birth. They established 
global existence and uniqueness of solutions using results from the theory 
of nonlinear evolution equations. Model (|1.1[ 1- (|1.4[) is a special case of the 
general model treated in [TJ, however, in [7J qualitative questions were not 
addressed. In contrast to [TJ, our paper focuses on the existence and local 
asymptotic stability of equilibrium solutions of system (|1.1[) - ()1.4|) with par- 
ticular regards to the effects of the distributed states at birth (previously 
we addressed simpler models without hierarchical structure in [12 | [M ) [15 ] ) . 
Earlier, in [TI5] Henson and Cushing studied continuous age-structured hier- 
archical models. They compared models with scramble versus contest com- 
petition for a limited resource. In particular the equilibrium levels for the 
two modes of competition were analyzed. Crucially however, the models in 
|19) incorporate vital rates that do not depend explicitly on the structuring 
variable. 

In Section [3J we will establish conditions for the existence of positive 
steady states of our model. The question of the existence of non-trivial 
steady states is difficult mainly for two reasons. Firstly, due to hierarchy in 
the population related to individual size, individual mortality and fertility 
depend on the size specific environment E. This environmental feedback 
yields an infinite dimensional nonlinearity in the model equations, in con- 
trast to Gurtin-MacCamy type models [IS] , where the vital rates depend 
on a weighted total population size, or on a finite number of such variables. 
Secondly, as individuals may be recruited into the population at all possible 
sizes, a recruitment operator of infinite rank arises. This means that the 
steady state equation cannot be solved explicitly. In [T2] we overcame this 
issue for a simpler model where the model ingredients only depended on the 
total population size by using results from the spectral theory of positive 
operators. Unfortunately this approach cannot be extended to the model 
considered here. Therefore we devise a different approach, based on fixed 
point results for nonlinear maps in conical shells of a Banach space, see 
[H [TU] . This approach was used before to treat age-structured models (also 
with diffusion), see |24M26) . where every individual enters the population at 
a single state, namely at age zero. The method is based on the construction 
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of an appropriate nonlinear map that requires the implicit solution of the 
steady state equation. However, the solution of the steady state equation of 
our model is not available. Therefore we need to construct an appropriate 
sequence of recruitment processes of finite rank for which we can solve the 
corresponding steady state problems. Then we show that the steady states 
constructed from the fixed points of the sequence of the nonlinear maps, 
have actually a convergent subsequence and we show that the limit point is 
actually a steady state of the original problem (|l.ip - (|1.4[) . 

In Section [3] we focus on the asymptotic behavior of solutions of the 
model. A positive quasicontraction semigroup describes the evolution of 
solutions of the system linearized at an equilibrium solution. We establish 
a regularity property of the governing linear semigroup in Proposition 13.71 
that allows in principle to address stability questions of positive equilibrium 
solutions of (|l.ll) - (|1.4[) . However, even the point spectrum of the linearized 
semigroup generator cannot be characterized explicitly via zeros of an as- 
sociated characteristic function. This is because the eigenvalue equation 
cannot be solved explicitly due to the infinite dimensional nonlinearity in 
the original model and the very general recruitment process. We will over- 
come this issue by using compact positive perturbations of the semigroup 
generator and rank one perturbations of the general recruitment term. This 
allows us to arrive at stability /instability conditions for the steady states of 
our model. 



2. Existence of positive equilibrium solutions. In this section we will 
discuss the existence of steady states of model (|1.1[) - (|1.4[) . We define the 
nonlinear operator 

by 

= 2£ (7( s : QMs)) + H(s, E{s 7 q))q(s) - J f3(s, y, E{y, q))q(y) dy, 

(2.6) 

where 

r-ra 

Q = K(s)q(s) ds. 
Jo 

It is clear that p* € T4^ 1,1 (0, m) is a steady state of ()l.ip - (|1.4j) if and only if 
= and p*(0) = 0. 
Our aim here is to apply a fixed point result, see e.g. Theorem 12.3 in 
4] or Theorem A in [21]. Its proof uses the Leray-Schauder degree theory 
for compact perturbations of the identity in infinite dimensional Banach 
spaces. 
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Theorem 2.1. Let (X, || • \\x) be a Banach space, K, C X a closed convex 
cone and K, r = K,C\B r {Q), where B r (0) denotes the ball of radius r centered 
at the origin. Let <f> : /C r — > K, continuous such that <S>(JC r ) is relatively 
compact. Assume that 

1. Qx^Xx for all \ \x\ \ x = r, A > 1. 

2. There exists a p e (0,r) and k € JC \ {0} such that 

x — <j>x ^ Xk for all \\x\\x = p, A > 0. 

Then <I> has at least one fixed point x* in the shell 

Sp.r = {x e K, : p < \\x\\ x < r}. 

2.1. The case of a separable fertility function. First we show in Theo- 
rem [52] that for a finite rank fertility process the problem admits a positive 
steady state under biologically meaningful conditions on the model ingredi- 
ents. Later we treat the general case in Theorem 12.31 

Theorem 2.2. Assume that 

i 

f3(s,y,E(y)) =J2PMPj(v,E(y)) (2.7) 
i=i 

with continuous functions (3, (3 and there exists a j G {1, 1} such that 

r-rn 

\ ^■( s ,0)F J -( S ,0,0)ds> 1, (2.8) 
Jo 

where is the zero function and 

'p(r,«W) + T.(r,P)J 



Lei F be a bounded and measurable function that satisfies 

F(s,H(s),P)>F j (s,H(s),P) 
for all j , s € [0, m], if 6 £+(0, m) and P > 0, and c be a constant such that 

i 

K(s)>cJ2f3 k {s,H{s)). (2.9) 

Suppose that there exists an R > suc/i £/ia£ /or a/Z (i?, P) G -^+(0, m) x R+ 
with ||P||li + P > R we have 

/ K(s)F(s,H(s),P)ds < c. (2.10) 
Jo 

Then the model (|l.l|l - (jl.4|l admits a positive steady state p*. 
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Proof. Let X = L'^m) © I 1 with norm || ■ \\ x = \\ ■ \\ L i + \\ ■ 
We will use the notation x=(H,P), where H G L^C^m) and 
P € I 1 or i= (H, P° ,P') where P° is just the first component of P and 
P' = (P 1 , P 2 , . . . ). We consider elements of l n to be in I 1 by the trivial 
embedding. We denote by K. = (i 1 (0, m) © the positive cone of X 
which is closed and convex and denote K r = K fl -B r (0), where r has yet to 
be chosen. Without loss of generality we may assume that the indices have 
been assigned such that condition (|2.8[) holds for j = 1. 

We note that for a fertility function j3 of the form (|2.7[) the non-trivial 
time independent solution of model (|l.lj) - (|1.4l) can be found (if it exists) as 

i 

p.{a)=Y t HF j {a,E.{a),F^), (2.11) 

3=1 

where 

Pl= p j (8,E.(a))p.(a)d8, je{l,...,l}, 
Jo 

E*(s) — a / w(r)p*(r) dr + / w(r)p*(r) dr, (2-12) 

PS = / «(*)!>.(*) d*. 
Jo 

We define a nonlinear map : /C r — >• /C by = ($' l7 $2)? where 
$[(H,P)(s) 




Although at this moment we consider only a single fertility function /3 
(namely the finite sum of separable functions from equation (|2.7I0 . we use 
the superscript notation in anticipation of a sequence of such maps that 
we shall employ in the proof of Theorem [231 below. Note that (H, P) > 
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and since 

Krtff.PlW^Pla-lM^t^fs),? ), 

3=1 

also ${(H, P) e VF 1,:L (0, m). Likewise $2(-ff, P) > and the nonzero entries 
of $2(-ff,P) are in Mv" 1 . Thus $ maps indeed into /C and has a relatively 
compact image. It can be seen that (i?*,P*) 6 /C r is a fixed point of the 
nonlinear map if and only if the function defined via equations (|2.11[) - 
(|2.12|) is a steady state of problem (|l.i p -(|1.4 |) . 
Assume now that for (H, P°, P') S /C r such that 

\\(H,P ,P')\\ x = \\H\\ Ll +P° + \\P% 1 =r 

and for some A > 1 we have 

$'(ff,P°,P') = X(H,P°,P'), 

that is 

<S>[(H,P°,P') = XH, <f> l 2 (H,P°,P')=\(P°,P'). (2.13) 
The second equation in (|2 . 13[) can be written as 

Xp o = y2 pj K ( s ) Fj (s,H( S ),P )d Sl (2.14) 
Apl =E P '7 ft(a(sF;(a(s),P°)d«, (2-15) 

.7 = 1 Jo 



X pl = J2 P j l ( s ,H(s))F j (s,H(s),P°)ds. (2.16) 

It follows from equations (|2.13[) - (|2.16[) that we may assume that P J ^ for 
.7=0,1,2,...,/. 

From equations (|2.15p - (|2,16p we obtain 

j j /• m 

A H P 'H' 1 =EE P " / ^(s,i?( S ))^( S ,P( S ),P°)d S . (2.17) 
fc=U=i Jo 
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Combining equations (|2.14[) and (|2.17l) wc obtain 

/ l t>irt l 



|P'II/1 =P° 



<- rm <■ 

i=l J ° k=l 

I p rn 

Y,Pi K^Fj^His), P°)ds 
„-_i Jo 



\ 3 = 1 

This, combined with inequality (|2.9[) implies that 

cUP'Ibx <P°. 



(2.18) 



However, by choosing r = \\H\\ L i + P° + ||P'|| ( i > \\H\\ L i + P° > R 
sufficiently large and using condition (|2.10p , we have from equation (|2.14[) 

j r-ra 

P" <VP J / K(s)F j (s,H(s),P°)ds 

p rn 

<||P'||,i / K{8)F(8,H(8),P°)ds<c\\P'\\ t i<P , 

Jo 

a contradiction. Thus condition (J} of Theorem 12. H is established. 

Let us now define k = (0, (1, . . . , 1, 0, . . . )) G K, \ {0} with I + 1 entries 
1 and assume that for some A > and p > we have for all (H, P) with 
\\(H,P)\\ x =p 

(H > P)-* , (H > P) = AA, 

that is 

H - <f[(H, P) = 0, P - ^(If, P) = (A, . . . , A). 
The latter equation can be written as 

(I-B(tf,P°)).P = (A,...,A), 
where the (I + 1) x (I + 1) matrix B has elements Bqq = and 



(2.19) 



P>, 



p l (s,H(s))F l (s,H(s),P )ds for i = l,...,l. 



It follows from condition (|2.8[) and the continuity of /?i and Fi that Sn > 1 
for all 1 1 (H, P°) 1 1 = p for some small enough value p > 0. This renders the 
left hand side of (|2 . 19[) negative and yields a contradiction. Condition ([2]) 
of Theorem 12.11 is established and the proof is now completed. □ 
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2.2. General fertility function. Every fertility function f3 of the required 
regularity C 1 (in all its arguments, see (|1.5[) ) can be written as a limit of 
partial sums of separable functions 

i 

(3 l (s, y, E(y)) = ]T P l k {s)i3 l k {y, E(y)), (2.20) 

k=l 

with 

lim ||/3(-, .,£(•)) -j8'(-, -,£(•)) | U=o ([0>m]2) =0 

for every .E 6 L (0, m). This can be achieved by partitioning the interval 
[0, m] into / subintervals by y k = k 7 -j- , fc = 0, . . . , I and setting 

P l k+1 (s)=P^y l k ,E{y h )) and p l k+1 (y,E(y)) = X b *_ ll3/ <](y) 

for fc = 1, . . . , Z where Xr„; , u j i denotes the indicator function of the reSpeC- 
L^fc — l ' "feJ 

tive interval. Then the partial sums converge in the suprcmum norm as the 
number of subintervals increases to infinity and there is a uniform bound of 
the derivatives in the s-direction 

¥k\\cH[o,m]) < M (2.21) 

for all k, uniformly for arbitrary E 6 L 1 . The latter is because of the 
boundedness of all derivatives of [3. 

Theorem 2.3. Assume that for fixed model ingredients fx, 7, /3 and k there 
exists a bounded and measurable function b that satisfies 

b(s) > j3(s,y,E) for every s, y £ [0,m] and £eL|(0,m) (2.22) 

and there exists a R > such that for \ \H\\i,i + P° > R we have 

/ K(s)F b (s,H(s),P°)ds<c, (2.23) 
Jo 

where 

/i(r,if(r)) + 7 s (r,P°) , 1 b(x) 



dx, 



Jo I Jx j{r,P u ) 

and c satisfies 

K(s)>cb(s), (2.24) 

/or every s €E [0, m] . Moreover assume that there exists a separable under- 
estimator of the fertility, 

< Px(a)fc(y,E(y)) < p(s,y,E(y)) (2.25) 

for every s 7 y G [0,m], E G L^_(0,m) swc/i that f3\ together with F\ satisfies 

p rn 

/ P 1 (s,0)F 1 (s,0,0)ds > 1. (2.26) 
Jo 
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Then model (|1.1[) - (|1.4[) admits a positive steady state p* € W ' (0,m). 

Proof. Wc begin by replacing (3(s,y,E(y)) by 

P(a,y,E{y))-^{aMy,E(y))>0 

from condition (|2.25|) and then decompose this remainder as indicated in 
equation (|2.20[) . We also note that the decomposition of the j3 l functions 
in (|2.20p can be achieved in the way that there is a common first term for 
every I, i.e. we may write /3i(s)j3i(y, E(y)) = p[(s)f3[(y, E(y)) for every I. 
This gives a sequence of finite rank approximations 

i 

P l (s,y,E(y)) = ^{s)Hy,E{y)) +J2ti( s Wl(y,E(y))- 

k=2 

Thus, by Theorem 12. 2[ every <£>' corresponding to a finite rank approximant 
P l has a fixed point P J, ) in the common shell 

Sp.r = {x £ K : p <\\x\\ x < r}. 

Indeed, due to the uniform bound from inequalities (|2.22[) and (|2.23[) . the 
outer radius r can be chosen uniformly. The common lower radius p can 
be guaranteed since all approximants /3 l begin with a common first term 
for which condition (|2.8p from Theorem 12.21 holds. By the additional gain 
in regularity due to we have that E\ € W 1 ' 1 (0,m). By straightforward 
calculations, 



7(s,P ' 

F"(s,El(s),Pl) 



~/(s,P?> 1 ) 7 2 (s,P*°'') 



— P,( s ,^( s ),Pt) 



ls ( S ,P^)(p(s,El( S ))+ ls ( S ,P^)) 



F 3 ( S ,El(s), P[) 



7(s,iV ) 

This shows that • ,£7* ( • )> p l) G W 2 ' x (0,m) for all I and j = 1, . . . , I, 
and moreover, that this family is uniformly bounded. For every I the fixed 
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point yields a steady state of the approximate problem by 

I 

i=i 

Since pi is a linear combination of elements in W 2,1 {Q,m), it is itself in 
W 2,1 (0, m). The uniform bound on the derivatives in (|2.21[) implies that 



p* 



li 1 



dpi 


+ 


d 2 p\ 


ds 


L 1 


ds 2 



< M for all / e N. 



This means that {.fl}^ is a bounded set in W 2,1 (Q, m) which is compactly 
embedded in W 1 ' 1 (Q,m) (see e.g. Theorem 6.2 in [3]). Therefore the se- 
quence pi has a convergent subsequence, again denoted by pi, with limit 
point in W 1 ' 1 (0, to). This p* is the natural candidate for a positive steady 
state of model (|l.ll) - (|1.4[) . Next we show that 

lim \ \V(pi)\\ =0, 

where ^> is the nonlinear operator defined in (|2.6[) . Let ty l denote the non- 
linear operator corresponding to the partial sum fertility function f3 l . Then 



I* (pi 



< 



< 



[p(s, y, El (y)) -(3 l (s, y, E\ (y))] pi (y) dy 
l(s,.,El(.))-p l (s,.,El(.))] pi(-)\\ds 
\\p(s,-,Ei(-))-p l (s,-,El(-))\\ LOO ds 



ds 



< I KHz 

<K\\p(-,-,Ei(.))-^(.,.,Ei(.))\\ LO 

for some positive constant K. Finally we have 



ll*(p* 



* ( lim pi J 




lim * (pi) 




\ l— >oc J 


L 1 


I— >OQ 


L 1 



= lim 



o, 



where the second equality follows from the continuity of the operator ^. 
Thus is the desired steady state. □ 

Remark 2.4 We note that conditions (|2.22[) - (|2.26[) are natural and biolog- 
ically relevant. They are similar to the ones obtained in [53] for the existence 
of a positive steady state of a nonlinear age-structured model. For our model 
the introduction of a net reproduction function (which will be an operator) 
will be somewhat cumbersome and biologically less straightforward, see the 
next section. However, it is still shown that conditions (|2.22l) - (|2.26|) require 
that the growth rate of the population is larger than one close to the zero 
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steady state while the growth rate of the population is small (definitely less 
than 1) for large population sizes. 

3. Asymptotic behavior. In this section we will investigate the asymp- 
totic behavior of solutions. Our approach is based on a formal linearization 
around a steady state solution and on a careful spectral analysis of the lin- 
earized operator. More precisely, our goal is to establish conditions which 
guarantee that the growth bound luq of the linearized semigroup is negative, 
respectively positive. We will show that the growth bound of the linearized 
semigroup can be completely characterized by the spectrum of its generator. 
However, as we noted before, the main difficulty is that the eigenvalues of 
the linearized semigroup generator and therefore the spectral bound cannot 
be characterized directly via eigenvalues, in the general case. 

We note that the Principle of Linearized Stability has so far only been 
established for semilinear models, see e.g. jTEJ 03J HZ] , but not for general 
quasilincar equations, such as the one we treat in this paper. Therefore, for 
the remainder of the paper, we make the additional assumption that the 
growth rate is separable 

7 ( Sj P)=7i(s)72(P). (3.27) 

This is plausible from the biological point of view, as the growth rate is 
modulated by the total weighted population, equally for individuals of all 
sizes. Then the quasilinear problem (|l.ip - (jl.4[) can be written in the form 

^=g(p)Ap + F(p), p(0)=po, (3-28) 



where 



Ap= (7ip) 



o 

and the recruitment and mortality terms in equation are incorporated 
in the nonlinear operator F. Grabosch and Heijmans in |17j introduced the 
transformation 

T P (t) = / g(p(s))ds 
Jo 

and defined 

q(r) =p(t p (r)), forr>0, 

where t p is the inverse function of t v . It is then verified that q satisfies the 
semilinear equation 

^ = Ag(r) + B(q(r)), q(0) = p , (3.29) 
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with the same initial value as in (|3.28|) . where the nonlinear operator B is 
defined via B(q) = F(q)/g(q). This requires that g is a continuous, strictly 
positive and bounded function, as is guaranteed by our assumptions (|1.5[) . 
In [T7] Grabosch and Heijmans showed that 

1. Solutions of problems (|3.28p and p.29[) are in one-to-one correspon- 
dence with each other [T7J Theorem 3.4]. 

2. Problems (f3T28|) and ([3^29)1 have the same equilibrium solutions, and 
an equilibrium solution of (|3 . 28[) is stable if and only if it is stable for 
(pT29| pll Theorem 5.1]. 

In what follows, we make the assumption of a separable growth rate (|3.27[) . 
but keep the notation j(s,P) to avoid cumbersome expressions. 

3.1. Linearization around steady states. Given a stationary (time in- 
dependent) solution of system (|1.1[) - (|1.4[) . we introduce the perturbation 
u = u(s,t) of p by making the ansatz p = u + p» and we substitute this 
into equations f|l . ![) - f)l .4[) . Then we are using Taylor series expansions of 
the vital rates of the following form 

f(x,E) = f(x,E*) + f E (x,E*)(E - E*) + "higher order terms", 

to arrive at 

u t (s, t) + ( ( 7 (s, P„) + 7 p(s, P*)U(t)) (u(s, t) + P *(s))) s 
= -(m(s,-EO>.P*)) + Ve(s,E(s,p*))E(s,u)) (u(s,t) +p*(s)) 

(/3(s,y,E(y,p*)) + (3 E (s,y, E(y,p„))E(y,u)) (u(y,t) + p*(y)) dy, 

(3.30) 

where we have defined 

pin 

U(t) = / k(s)u(s, t) ds. 



Next we omit the nonlinear terms in equation (|3.30|) to arrive at the lin- 
earized problem 

u t (s,t) +~f(s,P*)u s (s,t) + ~/ s (s,P*)u(s,t) +j Ps {s,P :t )p :t {s)U (t) 

+ 7p(s,PMs)U(t) 
= -[i(s, E(s,p*))u(s,t) - fj, E {s, E(s,p*))p*(s)E(s,u) 



i3(s,y,E(y,p*))u(y,t) + E {s,y,E(y,p*))p*(y)E(y,u) J dy, 

(3.31) 

with the boundary condition 

7(0, P*)it(0, t) = 0. (3.32) 
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Equations (|3.31|) - (|3.32[) are accompanied by the initial condition 

u(s,0)=u (s). (3.33) 

The linearized problem (|3.31|) - (|3.33|> is treated effectively in the framework 
of semigroup theory. To this end we cast the linearized system (|3.31[> - (|3.33[> 
in the form of an abstract Cauchy problem on the state space L 1 (0,m) as 
follows: 

^-u=(A + B + C + V + J 7 )u, u(s,0) =u (s), 
at 

where 

Au = - 7(-, P«) u s with Dom(i) = {u 6 W 1,l (0, m) \ u(s = 0) = 0} , 

(3.34) 

Bu = - (ls(; P*) + H(; E(;p*))) u, (3.35) 

p in 

Cu = - (jp s (-,P*)p*(-) +7p(->-P* y*(0) / K iy)u(y)dy 



= ~P*{-) K(y)u(y)dy, (3.36) 
Jo 

Vu = - fl E (; E(-,p*))p*(-)E(; U-), (3.37) 



Tu= \ P(-,y,E(y,p*))u(y)dy+ /?£;(-, y, E(y,p*))p*(y)E(y,u) dy. 
Jo Jo 

(3.38) 

Proposition 3.5. The operator A + B + C + V + J- generates a strongly 
continuous quasicontractive semigroup {T(t)} t>0 of bounded linear opera- 
tors on L 1 ^, m), which is positive if the operator C + T> + J- is positive. 

We omit the proof of the above proposition as it can be established 
following similar results in Q31 [HI US] • 

3.2. Stability results for positive equilibria. In this section we estab- 
lish linear stability /instability results for a general steady state p*. We will 
treat the extinction steady state in the next subsection in detail. When ad- 
dressing stability questions the main difficulty is that, in general, the point 
spectrum of the semigroup generator cannot be characterized via zeros of an 
associated characteristic function. We will overcome this problem by using 
positive perturbation arguments. 

Theorem 3.6. Assume that 

k(s) / exp^-/ — — dx (—( — WT dyds< ^ 1 - 

Jo { Jy l{x,P*) J l{y,P*) 

(3.39) 
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Furthermore assume that 

p*(s) < 0, s G [0, to] and 3e > such that p*(s) ^ for a.e. s G [0,e], 

(3.40) 

M»,£(*,p*))<0, M»»V»-E(y»P*))>0, 2/,sG[0,to]. (3.41) 

Then the steady state p* of model (|1.1|| - (|1.4|) is linearly unstable. 

Proof. First we note that A + B + C generates a positive, irreducible and 
eventually compact semigroup if conditions (|3.40p hold true. Eventual com- 
pactness is easily shown, see Proposition 13. 71 To establish irreducibility we 
consider the resolvent equation 

R(X,A + B + C)h = u, heX + . 

The solution is obtained as: 

u(s) - r exp | - r -y.(*,P.)+ti*,E( X ,p.)) + X 1 _h(yl_ 

TT f J [ s i a {x,P*) + p{x,E(x,p*)) + \ A \ p4y) 

- u l cxp rX dx /^y dy ' 

(3.42) 

where 



U = I k(s)u(s) ds 
Jo 

g g exp {- g 7,(x.P0+Mg.y d ,} _MhL_ dydg 

1 + f Q m «(a) £ exp { - /; 7.(x.P.H^y. P 0) + A d ,} _p^L d , ds ' 

(3.43) 

The second equality is obtained from (|3.42[) by multiplying by k and inte- 
grating from to to. Hence for A large enough U > and u ^s> follows 
from condition (|3.40p . The irreducibility and eventual compactness of the 
semigroup imply that the spectrum a(A + B + C) is not empty, see e.g. Th. 
3.7 in Sect. C-III in [5], hence the spectrum a(A + B + C) contains at least 
one eigenvalue. Next we find the solution of the eigenvalue equation 

(A + B + C)u = Xu 



as 



f f J s (x,P*) + p(x,E(x,p*)) + X \ p*{y) A 
u(s) = -U exp <^ - / — — dx \ — — — dy. 

(3.44) 
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Then wc multiply the solution (|3.44[) by k and integrate over [0, m] to obtain 

U = -U k{s) 
Jo 

Jo I Jy 1{X,P*) J l{V,P*) 

We note that, if U = then equation (|3.44[) shows that u(s) = 0, hence 
we have a non-trivial eigenvector if and only if U ^ and A satisfies the 
following characteristic equation 

1 = K(X) d ^ 

r r_ r ls (x,p*) + n(x,E(x, P *)) + x \ My) 

<0 P l Jy 1(X,P.) dX fl(y,P*) V 

(3.45) 



JO 



It is easily shown that 



lim K{\) = 0, 

A— ^+oo 



therefore it follows from the assumption i^(0) > 1 (|3.39[) . on the grounds of 
the Intermediate Value Theorem, that equation ()3.45j) has a positive (real) 
solution. Hence we have 

< s(A + B + C), 

where s(A + B + C) stands for the spectral bound of the operator A + B + C. 
The operators V and T are positive if conditions (|3.4ip hold true. We have 
for the spectral bound (see e.g. Corollary VI. 1.11 in pTj ) 

< s(A + B + C) < s(A + B + C + V + T). 

Since the growth bound of the semigroup is bounded below by the spectral 
bound of its generator, the proof is completed. □ 

We note that the instability conditions (|3.41[) imply that mortality is a 
non-increasing function of the environment and fertility is a non-decreasing 
function of the environment. 

Next we establish conditions which guarantee that the equilibrium solu- 
tion is linearly asymptotically stable. To this end we establish first that 
the spectrum of the semigroup generator A + B + C + V + T consists of ei- 
genvalues only and that the spectral mapping theorem holds true. Then the 
growth bound of the semigroup equals the spectral bound of its generator. 
These follow however from the following result (see e.g. Corollary IV. 3. 12 

in El]). 

Proposition 3.7. The semigroup {T(t)}t>o generated by the operator 
A + B + C + V + J- is eventually compact. 
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Proof. We only sketch the proof here since analogous results for simpler 
problems can be found in [JJJ Q3J [TSJ US]. Due to the zero flux bound- 
ary condition and the finite maximal size, the operator A + B generates a 
nilpotent semigroup. The biological interpretation is that in the absence of 
recruitment the population dies out independently of the initial condition. 
C is a bounded linear operator of rank one, hence it is compact. It only re- 
mains to establish that the bounded linear integral operators T> and T are 
compact. These however, can be deduced using the Frechet-Kolmogorov 
compactness criterion (see e.g. Chapter X in [28]) from the regularity as- 
sumptions we made on the model ingredients, see the proof of Theorem 12 
in [T2] for more details. □ 
The previous result guarantees that stability is determined by the lead- 
ing eigenvalue of the semigroup generator A + B + C + V + J 7 , unless the 
spectrum is empty. In that case one further needs to establish positivity of 
the semigroup which guarantees that the growth bound coincides with the 
spectral bound, which by definition equals minus infinity. However, as we 
noted before, the eigenvalue equation 



cannot be solved explicitly. This is due to the infinite dimensional nonlinc- 
arity in the original problem and to the very general recruitment term. 

Compact perturbations do not change the essential spectrum of the semi- 
group generator. Therefore our approach to establish stability is to find a 
positive compact perturbation of the generator for which we can charac- 
terize the point spectrum via zeros of an associated characteristic function. 
To this end, for a separable fertility function j3, we introduce the following 
operators on L 1 (0, m), 



(A + B + C + V + J 7 - XT) u = 0, 



u(s = 0) = 0, 




We formulate our main stability result in the following theorem. 



Theorem 3.8. Assume that there exists a function 



f3(s, y, E(y,p.)) = /3i(s)/3 2 (y, E(y,p*)) 



such that 



P{s,y,E(y,p*)) < @(s,y,E(y,p*)), s,y <E [0,m]. 



(3.46) 
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Furthermore assume that the following conditions hold true 

H E {s,E(s,p*))<0, p E (s,y,E(y,p*)))>0, s,ye[0,m], (3.47) 
7p s (s,P*>*(s) + 7p(>,P*K(s) < 0, s G [0,m], (3.48) 

and i/ie characteristic function K^{\) given by equation (|3.54|) correspond- 
ing to f3 and to the modified operators T> and T 2 does not have a zero with 
non-negative real part. Then the stationary solution p* is linearly asymp- 
totically stable. 

Proof. We obtain the solution of the eigenvalue equation 

(A + B + C + V + T 1 +T 2 -Xl)u = 0, (3.49) 

as 



(3.50) 



u(*)=^ / fo(s,y)A(y)dy 
Jo 

+ U 2 f ti(s,y)f2(y)dy + U 3 I f X (s,y)f 3 (y)dy 1 



where 

XJ\ = / k(s)u(s) ds, JJ% = / w(s)u(s)ds, 
Jo Jo 



U%= I (3 2 (s,E(s,p*))u(s)ds, 
Jo 

fo\s, y) = exp {- f y 7.(x,P0+^(x,p0) + A ^ 

, -jPs{y,P*)p*{y) - ip{y,P*)p'*{y) , Pi(v) 
my) = -, — ' h\v) ~ 



i{y,P*) i(y,P*) 

Mi/) 



Pi (v) J™ P2(x,E(x, p*))p* (x) dx — p.E(y, E(y,p*) )p* (y) 



j(y,P*) 

We multiply equation p.50|) by k, w and by f3 2 and integrate from to to, 
respectively to obtain 

E7i(l + 011(A)) + U 2 a 12 {\) + U 3 a 33 (X) = 0, (3.51) 
f/i0 2 i(A) + C/ 2 (a 22 (A) + 1) + U 3 a 23 {\) = 0, (3.52) 
C/io 3 i(A) + C/ 2 a 3 2(A) + C/ 3 (a 3 3(A) + 1) = 0, (3.53) 
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where for i = 1,2, 3, 

aii(\) = k(s) / A (s,y)/i(y)dyds, 
Jo Jo 

/•in />s 

oai(A) =/ w(s) / fo(s,y)fi(y)dyds, 
Jo Jo 

pm p s 

a 3i (X)= p 2 (s,E(s,p*)) ti(s,y)fi{y)dyds. 



If A € a (A + B + C+V+F1+F2) then the eigenvalue equation (|3.49[) admits 
a non-trivial solution it hence there exists a non-zero vector (£/i, C/ 2 , ?7 3 ) 
which solves equations (|3.51[) - (|3.53|) . To the contrary, if (U\, [/ 2 , U3) is a 
non-zero solution of equations p.51|) - p.53|) for some A G C then (|3.50p 
yields a non-trivial solution u. This is because the only scenario for u to 
vanish would yield 

-Ui [ f^{s,y)h{y)dy^U 2 I fi(s,y)f a (y)dy 
Jo Jo 

+ U 3 [ ti(s,y)f 3 (y)dy, 



for every s G [0,m]. This however, together with conditions (|3.47[) - (|3.48[) 
would yield U\ = U2 = U3 = 0, a contradiction. Thus A G C is an eigenvalue 
of A + B + C+V+/F1+/F2 if and only if it satisfies the characteristic equation 

/ l + ou(A) oia(A) 013(A) \ 
i^(A) d =det 021(A) l + a 22 (A) a 23 (A) =0. (3.54) 

\ a 3 i(A) a 32 (A) l + a 33 (A) / 

Next we observe that conditions (|3.47p - (|3.48[) guarantee that both C, V and 
/F are positive operators. Therefore we conclude that A + B + C + V + J 7 
is a generator of a positive semigroup. Moreover, the operators (V — T>), 

(J-"f — J 7 ! ) and {F2 — are au positive and bounded. We have 
s{A + B + C + V + F) = s (A + B + C + V + jf + 

< s ( A + B + C + V + jf + J* + V - V + t\ - 7f + T P 2 ' - J"| 

A + B + C + V + F^ + ffj < 0, 

and the proof is completed. □ 

Example 3.9 The crucial assumption of the previous theorem is that the 
characteristic function Kg does not have a root with non-negative real part. 
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Here we only present an example when this condition may be easily verified. 
In particular, let us assume that 

K = Cl, W = C 2 , /3 2 = C3 

for some constants ci,C2,C3 > 0. In this special case the characteristic 
equation K^(X) = takes the simple form 

ti(s,y)(c 1 f 1 {y) + c 2 f 2 (y) + c 3 f 3 (y))dyds. 

Therefore there exists a unique dominant real eigenvalue, which is negative, 
if 

expi- — — da 

I Jy 1{X,P*) 

( Pi{y) (c 3 + c 2 /" P 2 {s, E(s,p*))p*(s) ds) 

x ^ i(y,P*) 

c 2 fi E (y,E(y,p*)) + cup s (y,P^)p^(y) + c 1 ^ P {y,P if ) 

7{y,P*) 

< i. 



dy ds 



The latter condition may be easily verified for fixed model ingredients. 

3.3. The extinction equilibrium. In this subsection we establish a sim- 
ple criterion for the stability/instability of the trivial steady state = 0. 
The stability of the trivial steady state is important from the biological 
point of view, as it can answer the question for example if a species can be 
introduced successfully into (or can invade) a new habitat. In case of the 
trivial steady state the eigenvalue problem can be written as 

(A + B + Tp - XL)u = 0, u(s = 0) = 0, (3.55) 

where A 6 C, the operators A and B are defined via (|3.34p ~ (|3.35[) with 
p* = (and P* = 0) and 



Jo 

We recall that in case of simpler scramble competition models the so called 
net reproduction function played a crucial role in the stability analysis of 
equilibrium solutions, see e.g. [13] - In particular, we managed to relate our 
stability results to a biologically meaningful model ingredient, namely the 
net reproduction rate. In case of hierarchical contest competition models 
this is less straightforward as we have seen. Nevertheless at least for a 
separable fertility function 
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we may define a net reproduction functional 

R : L\(0,m) -> K 
of the standing population p as 
R( P ) 

ex P { - f ^P)^E [Xl p)) ^| My)Ms,E(8, P )) dy ^ 

Jo 1[S,P) Jo I Jy l(x,P) J 

where P is the weighted population according to (|1.4[) . The value of the 
functional R is the expected number of offspring to be produced by an 
individual in her lifetime. Individuals of size y are produced at a rate /3i 
and need to survive from size y to size s to reproduce. We also note that 
R{p*) = 1 is a necessary condition for the existence of a positive steady 
state p* of our model. It is however not a sufficient condition in contrast to 
scramble competition models, i.e. when density dependence is incorporated 
via finite dimensional nonlinearities in the vital rates. 

Proposition 3.10. Assume that there exists a function j3 such that 

p l {s,y)=P[{ S )P l 2 {y)<P{ S ,y,Q) 

and Rpi (0) > 1 (where Rpi stands for the net reproduction functional corre- 
sponding to the fertility function j3 l ). Then the trivial steady state is linearly 
unstable. On the other hand if there exists a function f3 w with 

and Rpu(0) < 1, then the trivial steady state is linearly stable. 

Proof. The eigenvalue equation (|3.55j) has a nontrivial solution if and only 
if A satisfies the characteristic equation 

1 = A-'(A) d ^ f 

r r 7s (,,f») + ^% Pt )) + A i p\{y) 

P l J y i(*,P*) )i(y,P*) y 

(3.56) 

It can be shown that lim K l (X) = 0. Therefore, if Rai (0) > 1 holds true, 

A— 5- + 00 

then the characteristic equation (|3.56[) admits a positive solution A. Since 
A + B + J- pi generates a positive semigroup, and the operator J-p — Tpi is 
positive we have < s(A + B + T m ) < s(A + B + J-p) and the instability 
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part follows. Similarly, if Rpu(Q) < 1 holds true, then the characteristic 
equation 1 = K U (X) does not have a solution with non- negative real part. 
Therefore on the grounds of Proposition 13.71 we have 

s(A + B + T fj + J>„ - Tp) < s(A + B + J» < 0, 

and the stability part follows. □ 

4. Concluding remarks. Hierarchical structured population models are 
important from the application point of view as they describe the evolution 
of populations in which individuals are experiencing size-specific environ- 
ment. This is the case for example in a tree population or in a cannibal- 
istic population. They pose a greater mathematical challenge though than 
scramble competition models. In the recent papers [14j [15] we started to 
investigate the asymptotic behavior of hierarchical structured partial diffe- 
rential equation models with one state at birth. Nevertheless, the question 
of the existence of positive equilibrium solutions for hierarchical models re- 
mained an open question up to our knowledge. In this work we formulated 
biologically relevant conditions for the existence of positive steady states of 
a very general model. 

In the second part we focused on the stability of equilibria. As we have 
also seen earlier in [21 [15], the main difficulty to establish easily verifi- 
able stability /instability conditions for positive steady states of hierarchical 
structured population models is that the point spectrum of the semigroup 
generator cannot be characterized explicitly via zeros of an associated char- 
acteristic function. In [14j [15] we devised a dissipativity calculation to show 
that the linearized semigroup has negative growth bound, and therefore 
the steady state is stable. The stability conditions we obtained in that way, 
though, were extremely restrictive. Therefore in this work we devised a new 
approach to establish stability by using results from the theory of compact 
and positive operators. 

The separability assumption made in (|3.27[) is only required for the result 
of Grabosch and Heijmans |17j to be applicable, and not to prove any of 
the results concerning the linearized semigroup T(t). A Principle of Lin- 
earized Stability for fully nonlinear quasilinear hyperbolic equations remains 
a problem for future investigations (see also (22] [25j [27] for related results) . 
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